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Abstract 

We explore the properties of dark energy from recent observational data, including the Gold Sne 
la, the baryonic acoustic oscillation peak from SDSS, the CMB shift parameter from WMAP3, the 
X-ray gas mass fraction in cluster and the Hubble parameter versus redshift. The AC DM model 
with curvature and two parameterized dark energy models are studied. For the AC DM model, we 
find that the flat universe is consistent with observations at the la confidence level and a closed 
universe is slightly favored by these data. For two parameterized dark energy models, with the 
prior given on the present matter density, 0,mo, with = 0.24, = 0.28 and ^ImO = 0.32, 
our result seems to suggest that the trend of ^ImO dependence for an evolving dark energy from a 
combination of the observational data sets is model-dependent. 
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I. INTRODUCTION 



The present cosmic accelerating expansion has been confirmed by various observations, 
including the Type la Supernovae (Sne la) [l, 13, S 0, 0, ^, CMB Q, S S Q, 111 and 
large scale structure (LSS) Q, Q, etc. In order to explain this observed phenomenon, it is 
usually assumed that there exists, in the universe, an exotic energy component with negative 
pressure, named dark energy (see [l4, 15, 16, 17 1 for recent reviews), which presumably 



began to dominate the evolution of the universe on 



dark energy is the cosmological constant A 



18 
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y re cently. The simplest candidate of 



2l| . It fits the observational data very 



well, but at the same time, it also encounters two problems, i.e., the cosmological constant 
problem (why is the inferred value of cosmological constant so tiny (120 orders of magnitude 
lower) compared to the typical vacuum energy values predicted by particle physics?) and the 
coincidence problem (why is its energy density comparable to the matter density right now?). 
Therefore, some dynamical scalar fields, such as quintessence 



22 
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quintom 



2J], phantom p5| and 



26( 1 ■ etc, are suggested as alternative candidates of dark energy. One of the features 



of these scalar field models is that their equations of state parameter, w, which embodies 
both gravitational and evolutionary properties of dark energy, is evolving with the cosmic 
expansion. 

On the other hand, the growing number of dark energy models has prompted people to 
adopt a complementary approach, which assumes an arbitrary parametrization for the equa- 
tion of state w{z) in a model-independent way and aims to reconstruct the properties of dark 
energy directly from observations. Currently, there are many model independent parameter- 



izations (see for example. 
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321]). In general, using these parameterizations 



and the observational data, one can determine the present value of w and whether it evolves 
as the universe expands, in particular, whether the phantom divide line (PDL) is crossed. 



In this regard, Nesseris and Perivolaropoulos 33| has used the Chevallier-Polarski-Linder 
parametrization w{z) = wo + Wiz/{1 + z) [28^ to explore the properties of dark energy with 
some observational data (including new Gold Sne la, SNLS Sne la, CMB, BAO, the clus- 
ter baryon gas mass fraction(CBF) and 2dF galaxy redshift survey ( 2 dFGRS) ) and found 
that the Gold data set mildly favors dynamically evolving dark energy with the crossing of 
the PDL while the SNLS does not, and the combination of CMB+BA0+CBF+2dFGRS 
mildly favors the crossing of PDL only for low values of VLmo ( ^mo ^ 0.25) prior consid- 
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ered and with a higher prior matter density the evolving features of dark energy becomes 



weaker and weaker. Similar trend of Qmo dependence was found using the model |34| . 
w(z) = i±£ di±M2(l±£) _ I ^ith the CMB and BAO. However, constraints from a combina- 



tion of the supernovae and other observational data has not been analyzed in Ref. 33|], and 



although that of the Sne and CMB+BAO was examined in Ref. 3J], but the marginalization 



was considered only for Q^o = 0.28 ±0.03 prior. Therefore, it remains interesting to see what 



happens to the conclusions reached in Refs. 



33, 



3J], when the combination of all observa- 



tional data is analyzed for different Qmo prior considered. The present paper aims to fill the 
gap. We discuss the constraints from the combination of different observational datasets. 
Besides the data sets of Sne la, BAO and CMB, in our analysis we add the datasets of the 
X-ray gas mass fraction in cluster and the Hubble parameter versus redshift. Firstly the 
ACDM model with curvature is discussed. Then, two parameterized dark energy models: 
w{z) = Wo + Wiz/{1 + z) and w{z) = ^-^^^^-^^^-^ ^ 1, are studied to see if the properties 
of dark energy thus reconstructed are model- independent. 

II. THE OBSERVATIONAL DATA 
A. The Gold Sne la data 

The Sne la data considered in this paper is the 182 Gold set. This set was released by 
Riess et al. {4]^ with a consistent and robust manner. It consists of 119 previously published 
data points [3], 16 points discovered recently by the Hubble Space Telescope (HST) and 



47 points from the first year release of the SNLS dataset [5[. For these Sne la, the data is 
released with the form of distance modulus fx, which is relative with the luminosity distance 
through 

fi{z,Ho,Pj) = 5\og^o[d^{z,pj)]+M . (1) 

Here for a flat universe = {1 + z) dz' / E{z' ,pj) with E'^{z) = H'^{z)/H^ = Qrno + 
(1 - n^o)exp[3 /o j^(l + w{z'))], M = M - 5logiQ{Ho), Hq = lOO/i, M is the absolute 
magnitude of the object and pj denote the model parameters of dark energy. The constraints 
on the dark energy models from Sne la data can be obtained by the maximum likelihood 
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method, so the best fit values for model parameters is determined by minimizing 



Xsne[J^O,Pj) - U [^) 



S,2 



B. The baryonic acoustic oscillation peak 

From the large scale correlation function of luminous red galaxy in the Sloan Digital Sky 
Survey (SDSS), Eisenstein et al. [3] found a baryonic acoustic oscillation peak, which is 
consistent with the prediction from the acoustic oscillation in the primordial baryon-photon 
plasma at the recombination. Thus this peak remarkably confirms the Big Bang cosmology. 
Meanwhile it also provides a ruler to constrain the dark energy models, which can be used 
usually by a dimensionless parameter A, 

'in dz 



2/3 



(3) 



E{zi,pjy/^[zi Jq E{z,pj) 
for a flat universe, where zi = 0.35 and A is measured to be A = 0.469±0.017. The parameter 
A is model-independent and clearly independent of the value of h too. By minimizing 

we can obtain the constraints from BAO. 



C. The CMB shift parameter 



For the CMB data, we use the shift parameter R to research the properties of dark energy, 



which can be expressed as 



35| 



RiPj) = V^mO 



dz 



(5) 



'0 E{z,pj) 

for a fiat universe, where Zr = 1048[1 + 0.00124((]fe/i2)-o.738] ^ gii^oh^^], gi = 
0.078(fib/i2)-o-238[i + 39.5(^^/^2)0.763]-! ^ 0.56[1 + 21.1(1^2)1.81]-! [ggj _ calculate 

Zr, we let Qbh^ = 0.024. The resuhs of three-year WMAP data n| give R = 1.70±0.03 Q- 



Let us note that the quantity R from CMB measurement is dependent on the value of VL^qK^, 
thus when using this parameter it is required to marginalize over Vl^fioh or assign some spe- 
cific value to h. In our discussion we give a prior value h = 0.72. We then place constraints 
on cosmological models using this shift parameter by minimizing 

,2 ^ _ [RiPj) - 1-70]' 



XcMsiPj) 



0.032 



(6) 
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D. The baryon gas mass fraction of galaxy cluster 



Under the basic assumption: the baryon gas mass fraction in cluster is constant, a compar- 
ison of the gas mass fraction of galaxy clusters, fgas_= Mgas/M^ot, can be used to constrain 



the cosmological models. Following Allen et al. [38|, l39| we fit the fgas data to a model 
described by 



fgas i^jPj) 



61^5(2/1)3/2 



3/2 



(7) 



;i + o.i9/iV2)a^o 

where 6 is a parameter motivated by gas dynamical simulations, d-A = d^/{^ + is the 
angular diameter distance, d^^^^ is the angular diameter distance corresponding to the 
standard cold dark matter (SCDM) universe (fimo = 1 for a fiat universe). Following 
Nesseris and Perivolaropoulos [33], we define A 



, „ — and treat it as a nuisance 

parameter. Using the theoretical method given in [33] to marginalize over A, we can obtain 
the constraints by minimizing 



whereA = E.(^)',i? = E 

J gas,i 

26 cluster data points given in Ref. 



Jga 

rmod robs 
J gas Jgas,i 



^ A ' 

/; 



(8) 



J'jas 

39| are used. 



^SCDi 



13/2 



gas 



E. The Hubble parameter data 



Based on differential ages of passively evolving galaxies determined from the Gemini Deep 



42] obtained 



Deep Survey [40|] and archival data [41j at redshift < 2 < 1.8, Simon et al. 
9 data points of H{z) at redshift Zi, which can be used to test the cosmological models by 
minimizing 

[Hobs{zi) - Hth{zi, HQ,Pj)f 



a 



(9) 



Hi 



Recently these 9 Hubble parameter data points have been studied extensively by many 



authors 



42 



431 ] ■ however it doest not provide a tight constraint on dark energy models. 



Thus the constraints on cosmological models from a combination of above discussed 
observational datasets can be obtained by minimizing 

X^{Ho, Pj ) = Xlne {Ho ,Pj)+ xIaO {Pj ) + xImB {Pj ) + Xjgas (Pj) + Xh{Ho, Pj) . (10) 
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Since we are interested in the model parameters, Hn becomes a nuisance parameter and is 
marginalized by a theoretical method given in Ref. [4J] in calculating Xsne Xh- 



III. RESULTS 



The AC DM model with curvature is firstly discussed with the observational data. The 
results are shown in Fig. 1 for a combination of the above discussed observational data 
sets. At a 95.4% confidence level we obtain fi^o = 0.29lo;o4 and fifco = -0.016lo:o29 with 

= 196.8. It is easy to see that a spatially flat universe is consistent with the observations 
at a 68% confldence level and a closed universe is somewhat favored by these data sets. 

TABLE L The best-fit data of Modi with prior value O^o- In the Table G+C+B, 
G+C+B+f and G+C+B+f+H represent the Gold+BAO+CMB, Gold+BAO+CMB+fgas and 
Gold+BAO+CMB+fgas+H(z) respectively. 



mO 



0.24 



Wo 



Wl 



XMin 



n 



mO 



0.28 



Wo 



Wl 



XMin 



mO 



0.32 



Wo 



Wl 



XMin 



Gold 
Gold+BAO 
CMB+BAO 
Gold+CMB 

G+C+B 
G+C+B+f 
G+C+B +f+Hl 



-1.28] 

-1.45 

-1.47 

-0.90 

-1.04 

-1.07 

-1.07^ 



1.59 
1.63 

52 

55 

85 



42 
37 

35 
31 

35 
1.30 

1.30 
1.30 



3.321:^ 
1.48j 



0.27 
0.63 
0.68 
0.71 



156.5 

il 158.1 

\il 0.001 

fs 161.8 

il 165.1 

fi 191.1 

■II 200.2 



-1.38 
-1.30 
-1.01 
-1.05 
-1.05: 

-i.io: 

-1.09" 



-0.67 
-0.72 

-0.53 
-0.59 

-0.94 
-0.71 

-0.44 
-0.40 

-0.39 
-0.34 

-0.40 
-0.32 

-0.36 
-0.32 



2.75 

2.43 
0.61 
0.81 
0.79: 
0.88: 
0.85 



-f3.52 
-3.73 

-f3.07 
-3.40 

4-1.40 
-6.08 

4-0.81 
-2.08 

-0.74 
-1.91 

-0.71 
-2.08 

-0.70 
-1.84 



156.5 
156.8 
0.001 
158.6 
158.6 
185.9 
194.9 



-1.49 



-fO. 



4-0. 



-1.11 



-o.6oto: 

-1.22l[]; 



-1.06 



4-0. 



-1.13 



+0. 



-1.12 



4-0. 



79 
85 

65 
64 

05 
71 

43 
56 

44 
36 

47 
35 

47 
35 



2.81 



4-4.13 
-4.62 



1.25 
-0.37 



4-3.50 



4.28 

1.82 
6.75 



-I 904-0.73 
J--^0-1.78 



0.96 



1.10 



1.06 



4-0.68 

-2.24 

4-0.62 
-2.54 

4-0.65 
-2.54 



156.6 
160.7 
0.001 
157.4 
160.^ 
189.6 
199.1 



Then we study, in the spatially flat case, the observational constraints on the following 
two-parameter models considered in Ref. 33| and 3J] respectively. 



Modi : w{z) = Wo + Wi 



(11) 



,0 (\ 1+^ ^1 + 2^2(1 + 2;) 
Mod2 ■ w[z) = — 1 , 



l-DE 



(12) 
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TABLE II: The best-fit data of Mod2 with prior value 17„o- In the Table G+C+B, 
G+C+B+f and G+C+B+f+H represent the Gold+BAO+CMB, Gold+BAO+GMB+fgas and 
Gold+BAO+CMB+fgas+H(z) respectively. 







^mO = 0.24 








^mO — 0.28 






f^mO — 0.32 






Ai 


A2 




Ai 




A2 




Ai 


A2 




Gold 


-3.6ll|i 




156.5 


_o oq+4.97 
"-'•'-'^-5.53 




156.5 


o 1 7+4.98 
"^•^'-5.54 


1 97+2.26 
-1.95 


156.5 


Gold+BAO 






158.3 


-3.01^ 


^4.48 
-5.09 


i.3otl:- 


156.7 




0.62^?:?! 


160.6 


CMB+BAO 




41+0.72 


0.001 


-0.211 


-2.90 
-2.58 


"•^^-0.79 


0.001 


1.08t|^? 


n 19+0.89 
^•^^-0.88 


0.001 


Gold+CMB 


0.6711:^,° 


-n 11+0-57 


161.6 


-0.33j 


-1.82 
-1.96 




159.1 


1 OO+2.02 


'-'•3^-0.69 


157.4 


G+C+B 




n 04+0-56 


166.1 


-0.35j 


-1.65 
-1.77 


'-'•^-'^-0.55 


159.2 


-o.7oli:9l 


^•■^^-0.62 


160.7 


G+C+B+f 


-0.051}:^^ 


-0.46 


192.1 


-0.53j 


-1.58 
-1.69 


'-'•^^-0.54 


186.3 


-o.9ili:S 


^•^'-'-0.60 


189.2 


G+C+B+f+H 




n in+o-52 


201.5 


-0.47j 


-1.54 
-1.60 


n 23+°-^'^ 

'-'•^"J-0.51 


195.3 


-0.73l|;62 


n QQ+0.63 
'-'••^"J-0.56 


198.8 



where fi/)^; = Ai(l + 2) + ^2(1 + 2)^ + 1 — fimo ~ ^1 ~ ^2 3J]. In order to find out whether 
the constraints from the observations are dependent on the choice of the value of fimo, we 
give three prior values of fi^o with Vtmo = 0.24, fi^o = 0.28 and fi^o = 0.32. Tables 1 
and 2 display constraints on model parameters with 95% confidence level. From the best fit 
values given in the Tables one can see that the BAO + CMB alone indicates that the evolv- 
ing features of dark energy become weaker and weaker with a higher prior matter density, 



which is the same as obtained in Refs. 



33 



34| . however the trend is reversed for the Gold 



+ CMB, Gold+CMB+BAO, Gold+CMB+BAO+fgas and Gold+CMB+BAO+fgas+H(z). 
Meanwhile we find the effects of adding the f-gas and Hubble parameter data sets are 
not very significant, suggesting that the model parameters are strongly constrained by 
CMB+BAO+Gold. In Fig. 2 we show the constraints on model parameters from all the 
data sets. The upper and down panels show the constraints on Modi and Mod2 respec- 
tively. In Fig. 2 the red dot denotes the fiat ACDM model. This figure clearly shows that 
ACDM model is consistent with the observations at the 95% confidence level. In addition 
one can find that the best fit of model parameters is closest to the ACDM when Qmo = 0.24, 
which can also be seen from the Tables. 

In Fig. 3 we plot the evolutionary curves of w{z) for these dark energy models with 
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different prior values over fimo- The upper and down panels show the results of Modi and 
Mod2 respectively. In Fig. 3 the solid lines show the evolution of w{z) with the best fit 
values for the model parameters, and the dotted lines are for la and 2a error bars. For 
Modi the best fit curves show that the combination of the data sets considered in this paper 
favors an evolving dark energy and a crossing of the phantom divide line in the near past, 
and suggests that the present value of w is very likely less than —1. Remarkably, these 
conclusions are almost insensitive to the chosen value of matter density, in a clear contrast 
to those obtained in 33|] where it was found that, when the Sne data are not combined 



with other observational ones, the evolving property becomes weaker and weaker with higher 
matter density and the phantom divide may not be crossed with the increasing of matter 
density. With this being said, it should be pointed out, however, that higher values of Qmo 
lead to larger errors, especially at the 2a level, as can be seen from the Tables. The results 
for Mod2 are shown in the down panels of Fig. 3. Unlike Modi, the best fit curves show that 
the properties reconstructed depend on the chosen value of matter density. For Qmo = 0.24, 
a very mildly evolving dark energy is obtained, but with the increase of Qmo prior considered. 



the evolving feature o: 
to that found in Ref. 



dark energy becomes more evident. This trend is just the opposite 



34| for just BAO+CMB data. These discrepancies can also be found 



in Tables 1 and 2. However at the 2a confidence level for Modi and Mod2 the cosmological 
constant cannot be ruled out. In addition, we also find that the stringent constraint on w{z) 
happens around redshift z ~ 0.5, which is consistent with that obtained in Refs. 3l|, l45[ |. 

Finally, the evolution of the decelerating parameter q with the redshift is studied for these 
two parameterized dark energy models with Qmo = 0.28. The results are shown in Fig. 4 
with the left and right panel corresponding to the results of Modi and Mod2 respectively. 
In this figure the dashed lines are the results of ACDM model with Qmo = 0.28, the solid 
lines show the evolution of q{z) when Q^o = 0.28, the model parameters are the best fit 
values, and the dotted lines are for la errors. The figure shows that the present value of 
the deceleration parameter go is less than zero, indicating that the cosmic is undergoing an 
accelerating expansion, but for different models the value of go is different. 
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IV. CONCLUSION 



In this paper, we have reconstructed the properties of dark energy from recent observa- 
tional data, including the Gold Sne la, the baryonic acoustic oscillation peak from SDSS, 
the CMB shift parameter, the X-ray gas mass fraction in clusters and the Hubble parameter 
data. The AC DM model with curvature and two parameterized dark energy models are 
discussed. We find that a spatially flat universe is allowed by these data sets at the 68% 
confldence level, and a closed universe is slightly favored by the observations. For two pa- 
rameterized dark energy models, we give the priors on with Q^no = 0.24, Q^no = 0.28 and 
^mo = 0.32. For the spatially flat case, the constraints on model parameters and the evolu- 
tions of w{z) and q{z) are studied. The Gold + CMB +BAO give the strong constraints on 
model parameters. If Modi parametrization is used, the best flt curves in Fig. 3 show that 
the combination of the data sets considered in this paper favors an evolving dark energy, 
a crossing of the phantom divide line in the near past, and the present value of w being 
very likely less than —1. Remarkably, these conclusions are almost insensitive to the chosen 



value of matter density, in a sharp contrast to those obtained in [33| where the Sne data are 
not combined with other observational ones. However, the best flt curves in Fig. 3 indicate 
that the properties of dark energy reconstructed using Mod2 parametrization depend on 
the chosen value of matter density. For = 0.24, a very mildly evolving dark energy is 
obtained, but with the increase of Qmo prior considered, the evolving feature of dark energy 
becomes more evident. This trend is just the opposite to that found in Ref. [3^ for just 
BAO+CMB data. Therefore, our result seems to suggest that the trend of Qmo dependence 
for an evolving dark energy is model-dependent. It should be noted, however, that at the 
2a confidence level, the cosmological constant are allowed for both Modi and Mod2. 
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FIG. 1: The la, 2a and Sa confidence contours for a AC DM universe with curvature from the 
combination of Sne la, BAO, CMB, the X-ray gas mass fraction in clusters and Hubble parameter 
data. 
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FIG. 2: The constraints on two parameterized dark energy models from the combination of the 
Gold Sne la, BAO, CMB, the X-ray gas mass fraction in clusters and the Hubble parameter data. 
The upper and down panels, respectively, show the results of w^z) = wo + wiz/{l + 2;)(Modl) and 
w{z) = i±£;il±HMl±£) _ 1 (Mod2) with three different priors over O^q: ^mo = 0.24, = 0.28 
and ^ImO = 0.32. The red dot represents the ACDM model. 
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FIG. 3: The evolution of w{z). The upper and down panels, respectively, show the results of Modi 
and Mod2 with Clmo = 0.24, Cl^no = 0.28 and ^^^0 = 0.32. The solid line shows the evolution of 
w{z) with the model parameters at the best fit values, and the dotted lines are for the la and 2a 
errors. 




FIG. 4: The behavior of q{z). The left and right panel, respectively, shows the results of Modi 
and Mod2 with ^1^0 = 0.28. The dashed line represents the ACDM with = 0.28, the solid 
line shows the evolution of q{z) with QmO = 0.28 and the model parameters at the best fit values, 
and the dotted lines are for the la error. 
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